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MiniBooNE, the Mini Booster Neutrino Experiment at Fermilab, will search for the appearance of electron 
neutrinos in a beam of muon neutrinos, and in the process will accumulate more neutrino interactions in the 
— 2 GeV energy range than all past experiments. This work presents preliminary results from the neutral 
current 7r° analysis at MiniBooNE, with comparisons of the NUANCE, NEUGEN, and NEUT neutrino cross 
section Monte Carlos. Also presented are comparisons of these three Monte Carlo simulations with MiniBooNE 
data taken during the running period from September, 2002, to September, 2003. 



1. INTRODUCTION 

MiniBooNE, the Mini Booster Neutrino Exper- 
iment at Fermilab, has been running smoothly 
and taking data for over a year now. The ex- 
periment intends to make a definitive statement 
about the neutrino oscillation signal seen by the 
Liquid Scintillator Neutrino Detector (LSND) 
Experiment at Los Alamos National Laboratory 
PP. Since both solar and atmospheric neutrino 
oscillations have been well established in recent 
years, confirmation of the existence of another os- 
cillation signal with its own distinct mass-squared 
splitting would indicate a need for further exten- 
sions to the Standard Model. 

On the way to answering the remaining ques- 
tion of the LSND signal in the neutrino oscil- 
lation landscape, MiniBooNE will collect more 
than 1 million neutrino interactions on pure min- 
eral oil (CH2). As there is very little experi- 
mental input to neutrino interaction cross sec- 
tions on heavy targets at MiniBooNE energies 
(0 < E„ < 2 GeV), the experiment's contribu- 
tions will be an important addition. 

Simulations of neutrino cross sections in the 
— 2 GeV region are fairly complex, and are mod- 
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eled differently by the three cross section Monte 
Carlos discussed in this work. The MiniBooNE 
cross sections group has worked in collaboration 
with the authors of the NUANCE 0, NEUGEN 
0, and NEUT gj Monte Carlos to provide model 
comparisons of neutral current 7r° (NC 7r°) pro- 
duction, comparing the three Monte Carlo simu- 
lations with each other and with the MiniBooNE 
NC 7T° data set. This is the first time kinematic 
comparisons of three different neutrino event gen- 
erators have been made to low energy neutral cur- 
rent 7T° production data. 

2. MiniBooNE OVERVIEW 

The Fermilab Booster accelerates protons to 
8 GeV kinetic energy in a 15 Hz cycle; 2 to 3 Hz 
are sent to MiniBooNE. This proton beam is di- 
rected into a 71 cm long beryllium target located 
inside a horn with a toroidal magnetic field. The 
horn focuses positively charged particles into a 50 
meter decay region. Any charged particles that 
do not decay in the 50 meters are stopped by an 
absorber located at the downstream end of the de- 
cay region. The neutrinos produced in the parti- 
cle decays then travel through approximately 500 
meters of dirt before reaching the detector. The 
beam at the detector, assuming no oscillations, is 
very pure - approximately 99.5% muon neutrinos 
with a mean energy of approximately 700 MeV 
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The detector consists of a sphere within a 
sphere, where the inner signal region is optically 
isolated from the outer veto shell [S]. The entire 
detector is instrumented with 8-inch Hamamatsu 
photomultiplier tubes (PMTs): 1280 in the signal 
region of the detector (10% coverage) and 240 in 
the veto region |7j. The entire detector is filled 
with approximately 800 tons of pure mineral oil 

E 

Beam arrives at the target in 1.6 /is bursts. The 
data acquisition system (DAQ) receives a signal 
from the Booster indicating that beam will be ar- 
riving imminently. This signal triggers the DAQ 
to open a window and record all detector activity 
for 19.2 /zs, starting 4.8 [is before beam arrives. 
Fig- U\ shows the distribution of events during the 
beam window. The top panel shows the 1.6 fis 
wide peak of the beam particles arriving during 
the beam window. 

With the application of only two very simple 
cuts, the background on the distribution can be 
drastically reduced. Requiring no more than 6 
hits in the veto region removes events with a co- 
incident cosmic muon occurring in the beam win- 
dow, shown in the center panel of Fig. ^ The 
further requirement of at least 200 hits in the 
main tank region, shown in the bottom panel, 
largely removes the chance of a muon decaying to 
an electron during the beam window. With only 
these two cuts, the signal to background ratio is 
greater than 1000:1. These cuts provide a very 
clean sample of beam-induced neutrino interac- 
tions to which further cuts may be applied for 
selection of particular interaction types. 
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3.1. External Predictions 

There are very few existing measurements of 
neutral current ir° production rates. A re-analysis 
of 1970's bubble chamber data from Gargamelle 
|9I10| is shown with cross section predictions from 
two neutrino cross section Monte Carlos in Fig. 
The predictions agree with the datum at its en- 
ergy, but there is currently no confirmation of 
agreement at other energies. Additional cross 
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Figure 1. The top panel shows recorded tank ac- 
tivity during the beam window with no cuts. The 
center and bottom panels respectively add a sim- 
ple veto cut and greater than 200 hits in the main 
volume of the detector. 



section measurements are available as ratios of 
neutral current to charged current 7r° production 
|9I11I12I13I14| . but these measurements can differ 
by up to factors of 3. 

The three neutrino cross section Monte Carlo 
(MC) simulations compared in this paper - NU- 
ANCE version 3, NEUGEN, and NEUT - share 
common theoretical inputs. All use the Llewellyn 
Smith free nucleon quasi-elastic cross section 
|15j . the Rein and Sehgal resonance cross section 
model and the standard deep inelastic scat- 
tering (DIS) formula for large momentum transfer 
|17|. There are nontrivial differences among the 
three though, viz., implementation of the Fermi 
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Figure 2. Neutral current single 7T° production 
cross section on protons as a function of neutrino 
energy. 



and cross sections to describe coherent it pro- 
duction processes [21] ■ NEUT rescales the Rein 
and Sehgal cross section down to the level pre- 
dicted by Marteau [23 , and of the three event 
generators, NEUGEN predicts the lowest coher- 
ent cross section. The coherent 7r° cross section as 
a function of neutrino energy is shown in Fig.|3for 
several models, including NUANCE, NEUGEN, 
and the Marteau model used to rescale NEUT. 
It should be noted that the NUANCE predic- 
tion shown here does not include pion absorption, 
which would lower the cross section prediction; 
absorption effects are included at a later stage in 
the event generation code. To demonstrate the 
wide variation in predictions for NC coherent 7T° 
production, the Paschos |2U and Kelkar/Oset [21] 
models are also shown. The two data points are 
from experiments that measured the cross section 
for coherent ir° production at 2 GeV and 3.5 GeV 
respectively. 



gas model for quasi-elastic interactions, joining of 
the resonance and DIS regions, and treatment of 
nuclear effects. The details of these MC simula- 
tions may be found in Refs. |2I3I4| . and also in 
other papers contributed to this set of conference 
proceedings |18ll9l20j . 

Neutral current 7r°'s are created by two main 
mechanisms: resonant production and coher- 
ent production. Resonant production is dom- 
inant, and occurs when a baryon resonance is 
excited and subsequently decays to a nucleon 
and one or more mesons, such as 7r°'s. For NC 
7T° resonance production, all three event genera- 
tors use a Rein and Sehgal-based model. Both 
NUANCE and NEUT use an axial-vector mass, 
m A = 1.1 GcV/c 2 , while NEUGEN uses m A = 
1.032 GeV/c 2 . 

In coherent NC tt production, a neutrino scat- 
ters from an entire nucleus rather than the indi- 
vidual constituents of the nucleus. In this pro- 
cess, very little of the momentum is transferred 
to the recoil nucleus, and the lepton and meson 
are produced in the forward direction. All three 
event generators use Rein and Sehgal kinematics 
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Figure 3. Low energy coherent neutral current 
single 7T° production as a function of neutrino en- 
ergy. A recent correction to the NUANCE gener- 
ator (version 3.004) results in a predicted coher- 
ent 7T° production cross section essentially identi- 
cal to that shown for NEUGEN. 
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3.2. Analysis Chain 

All Monte Carlo data sets used in this paper 
were produced by the same method. The Mini- 
BooNE neutrino flux was fed into each of the 
three cross section Monte Carlo simulations, and 
the generated events generated were then sent 
into the full MiniBooNE detector simulation. 
Events in each sample were reconstructed using 
the same algorithm. The NC ir° events were se- 
lected by applying the following cuts to each sam- 
ple: the background-reducing cuts discussed in 
section a fiducial volume cut (reconstructed 
event vertex R < 500 cm from detector center), 
and a required minimum energy for each of the 
gammas from 7r° decay (E 7 > 40 MeV). Using 
this standard analysis chain for each data set en- 
ables us to make a valid comparison of the three 
MC samples with each other and with data. 

Events which pass the event selection cuts 
have been reconstructed using a ir° fitting al- 
gorithm that fits events assuming there are two 
electron-like Cerenkov rings. The event kinemat- 
ics are determined by the direction and amount 
of Cerenkov light produced by the gammas from 
7T° decay. 

3.3. Monte Carlo Comparisons 

The momentum distributions for resonant and 
coherent 7r°'s which pass the NC w° event selec- 
tion cuts are shown for NUANCE, NEUGEN, and 
NEUT in the top panel of Fig. 0] normalized to 
the same number of events before cuts. The over- 
all rate for NC ir° production is very similar for 
NUANCE and NEUGEN, but quite a bit lower 
for NEUT. We also see that the NEUGEN mo- 
mentum spectrum is softer than that of both NU- 
ANCE and NEUT. The shape difference is more 
drastic when the distribution is separated into its 
resonant and coherent components, shown in the 
bottom left and bottom right panels of Fig.0J re- 
spectively. NEUGEN produces slightly more 7r°'s 
by resonance excitation than the other two Monte 
Carlos, and has a harder spectrum for coherent 
production. 

Fig.|S]shows the angular distribution of 7r°'s rel- 
ative to the beam direction for the three MC sam- 
ples. This distribution is sensitive to the produc- 
tion mechanism - coherent 7r°'s are much more 
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Figure 4. p^o distribution for Monte Carlo events 
passing MiniBooNE NC ir° event selection cuts: 
all NC 7T° (top), only resonant NC ir° (bottom 
left), only coherent NC ir° (bottom right). 



strongly forward-peaked than resonant 7r°'s, as 
shown in the bottom panels of this figure. NU- 
ANCE produces significantly more coherent 7r°'s 
than either NEUGEN or NEUT (bottom right 
panel). The lack of agreement among the three 
is not surprising, considering the many coherent 
7T° production models that exist and the dearth 
of data to rule them out. 

3.4. Data/Monte Carlo Comparisons 

Understanding the rate and kinematics of ir° 
production in neutral current interactions is criti- 
cal for MiniBooNE, since these events are a major 
background for the i/„ — > v e appearance search. 
The v e appearance search is a blind analysis; the 
box will not be opened until all the data have 
been collected. In the meantime, NC ir° events 
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Figure 5. cos 0^0 distribution for Monte Carlo 
events passing MiniBooNE NC 7r° event selection 
cuts: all NC 7r° (top), only resonant NC 7r° (bot- 
tom left), only coherent NC 7r° (bottom right). 
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Figure 6. Reconstructed invariant mass of beam 
triggers (open circles with statistical error bars). 
Fitted shapes are NUANCE-based parameter- 
izations of the contribution from background 
(dashed curve) and signal (NC resonant and co- 
herent 7r°'s, solid minus dashed curve). Errors 
are statistical only, and do not include system- 
atic errors associated with reconstruction. Note 
that some of the background events contain 7r°'s, 
so the peak near m^o in the dashed curve is ex- 
pected. 



can be used to test Monte Carlo predictions of 
event kinematics, and will also be used to mea- 
sure the neutrino neutral current 7r° cross section. 
A preliminary analysis of the mass spectrum and 
kinematics of 7r°'s produced in neutral current in- 
teractions observed in MiniBooNE is presented 
here with comparisons to the three Monte Carlo 
samples discussed above. 

Fig. shows the reconstructed invariant mass 
of beam triggers which satisfy the NC tt° event se- 
lection criteria. Cosmic-ray-induced activity dur- 
ing the beam trigger is eliminated by the method 
discussed in section [S] Further cuts on fiducial 
volume and energy of 7's from ir° decay are ap- 
plied to ensure good reconstruction of the event. 
The fitted curves are NUANCE-based parameter- 
izations, where the signal contribution arises from 



NC resonant and coherent single ir° events, and 
the background contribution comes from all other 
events. The background shape is determined by 
NUANCE version 3 simulations. A background 
peak near m^o is expected, since the background 
does contain some 7r°'s produced in final state in- 
teractions and some 7r°'s from multi-pion events. 
For approximately 1 x 10 20 protons on target, we 
extract 2509 ±113 signal- like NC tt° events in the 
fit to data in the preliminary analysis. The fitted 
mass peak agrees well with the nominal n° mass. 

The number of neutral current 7r°'s seen in data 
is also extracted in bins of calculated variables of 
interest (p n o, cos 9cm, and cos^o) to produce 
a distribution of each variable with the binned 
yields. The signal fraction in each bin is extracted 
via a fit to the invariant mass plot for events 
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in that bin. Unit area normalized distributions 
for each variable are compared with expectations 
from NUANCE version 3, NEUGEN, and NEUT. 
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Figure 7. Extracted ir° yields in p^o bins com- 
paring Monte Carlo and data, each normalized to 
unit area. Error bars on data are statistical only. 



The signal ir° momentum distribution for data 
is compared to NUANCE, NEUGEN, and NEUT 
in Fig. [7| Error bars shown on the data are sta- 
tistical only. Shape agreement between data and 
each of the three Monte Carlo samples is shown to 
be quite good. The fraction of 7r°'s reconstructed 
at high momentum drops off mainly due to the 
fall-off in the neutrino flux spectrum. 

Fig. [S] shows the cosine of the ir° center of mass 
angle - the angle between the direction of the 7r° 
in the laboratory frame and the decay axis of the 
7's in the center of mass frame. This quantity is 
related to the energy asymmetry of the 7's from 
7r° decay by cos 9cm — \ e!+e^ • Because the 7r° 
is pseudoscalar, it has no preferred decay direc- 
tion and the center of mass angular distribution 
should be flat. The spectrum is distorted, how- 
ever, for cos#cm >^ 0.7 because of the minimum 
energy requirement for each of the 7's in the de- 
cay; events where one or both of the 7's have 
less than 40 MeV are cut. All three Monte Carlo 
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Figure 8. Extracted ir° yields in cos 6cm bins 
comparing Monte Carlo and data, each normal- 
ized to unit area. Error bars on data are statisti- 
cal only. 



samples model the data well. This is especially 
important in the region where the spectrum falls 
off because these are the class of events that will 
mimic an oscillation signal. 

The angular distribution of 7r°'s relative to the 
beam direction is shown in Fig. As stated in 
section I5~3l this distribution is sensitive to the w° 
production mechanism; it may help us to deter- 
mine how much coherent production contributes 
to the overall rate of NC ir° production. The 
data appear to have a much lower contribution 
from coherent 7r°'s than all three Monte Carlos 
predict, as indicated by the deficit in the forward- 
most bin. 

4. CONCLUSIONS 

The work presented here is the first time sev- 
eral low energy Monte Carlo event generators in 
use by various experiments have been uniformly 
compared to experimental distributions of NC 7r° 
events. These preliminary comparisons to the 
MiniBooNE data are both interesting and encour- 
aging. The shapes of kinematic distributions are 
surprisingly similar despite their differing theo- 
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Figure 9. Extracted ir° yields in cos 9 v o bins com- 
paring Monte Carlo and data, each normalized to 
unit area. Error bars on data are statistical only. 



retical inputs. The exception is in the ir° angular 
distribution relative to the beam direction, where 
there is an expected difference due to the different 
coherent contributions among the three Monte 
Carlo simulations. The deficit in the forward- 
most bin of Fig.[5]seems to indicate a much lower 
contribution from coherent 7r° production than 
any of the MC simulations predict; more data 
are needed, however, before any conclusions may 
be drawn. There is a serious lack of NC 7r° cross 
section data at low energies for both resonant and 
coherent NC 7r° production; MiniBooNE's com- 
plete collection of NC tt° data, which will ulti- 
mately be ~ 10 times the data presented here, 
will provide an important measure of this cross 
section in the — 2 GeV region. 
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